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We investigate the properties of spontaneous currents generated at surfaces and interfaces of d-wave 
superconductors using the self-consistent quasiclassical Eilenberger equations. The influence of the 
roughness and reflectivity of the boundaries on the spontaneous current are studied. We show 
that these have very different effects at the surfaces compared to the interfaces, which reflects the 
different nature of the time reversal symmetry breaking states in these two systems. We find a 
signature of the "anomalous proximity effect" at rough d-wave interfaces. We also show that the 
existence of a subdominant order parameter, which is necessary for time reversal symmetry breaking 
at the surface, suppresses the spontaneous current generation due to proximity effect at the interface 
between two superconductors. We associate orbital moments to the spontaneous currents to explain 
the "superscreening" effect, which seems to be present at all ideal d-wave surfaces and interfaces, 
where dxy is the favorite subdominant symmetry. 



I. INTRODUCTION 

The qualitative difference between d-wave and conven- 
tional (s-wave) superconductors is the intrinsic 7r-shift of 
the order parameter phase between the crystallographic 
a and b directions The most remarkable effects due 
to it are reahzed near surfaces and interfaces of such su- 
perconductors. In particular, depending on the direction 
of the crystalline axes on both sides of a Josephson junc- 
tion, the equilibrium can be achieved not only at phase 
difference (f) — 0, like in conventional Josephson junc- 
tions, but ai (j) — TT as well ("TT-junctions") ||,||. There- 
fore, a frustrated ring (a ring with an odd number of 
TT-junctions) will have a time reversal symmetry break- 
ing ground state, which supports a spontaneous mag- 
netic flux ±<i>o/2, where = ir/e is the flux quantum 
(throughout this article we take h = ks = Experi- 
ments [H-|5] confirmed this prediction, thus establishing 
dj.2_y2 pairing symmetry of high- Tc cuprates. 

The TT-junctions are not specific to unconventional or- 
bital symmetry (being first predicted |^ and recently re- 
alized Q in SFS structures). They do not either provide 
the only way that T-breaking states can appear in d- 
wave superconductors. We can distinguish between two 
situations: a doubly connected geometry, like a ring with 
TT-junctions, and a simply connected one, like a single 
boundary between a d-wave superconductor and another 
differently oriented d-wave superconductor, s-wave su- 
perconductor, or vacuum. In the former case, the spon- 
taneous flux is quantized (in units of $o/2) 0. There- 
fore, the spontaneous currents are always present and 
flow through the whole ring (within the screening length 
from the surface which is, in high-Tj, cuprates, of order 
1500 A). In the latter case, on the other hand, there is 
no quantization condition for the spontaneous flux and it 
can take arbitrary values . The spontaneous currents 



can be absent or confined to a much narrower area near 
the surface/interface itself, as we will see later. 

A natural description of T-breaking states near the 
surfaces and interfaces of d-wave superconductors can be 
given in the language of Andreev bound states |^ . They 
are formed by off-diagonal scattering of quasiparticles by 
a spatially inhomogeneous pairing potential, A(r). Off- 
diagonal means that the reflected quasiparticle changes 
the branch of the dispersion law (particle to hole and vice 
versa), so that electric charge 2e is being transferred to 
or from the condensate. Remarkably, Andreev reflection 
conserves spin (exactly) and momentum (with accuracy 
E/Ep, E being the quasiparticle energy measured from 
the Fermi level, Ep)- As a result, the reflected hole is 
sent back with almost the same group velocity as the 
incident electron, and will therefore retrace its path (up 
to distance ~ — vp / E) in the clean limit (large elastic 
scattering length). 

Mathematically, the Schrodinger equation for the 
(quasi-) electron wave function is now replaced by a 
matrix Bogoliubov-de Gennes equation |10[| for two- 
component quasiparticle (bogolon) wave function ^'(r) = 

("(rj)- order parameter A(r) and its complex 

conjugate play the role of the off-diagonal components 
of the matrix scattering potential. If we neglect the 
self-consistency condition, which expresses A(r) through 
w(r), ^(r), the equations are easily solved, and the po- 
sitions of Andreev levels, E, are found from the Bohr- 
Sommerfeld quantization condition for a quasiclassical 
traj ectory bounded by Andreev reflections at points L, R 



/ p^,^SE)dl- p,^jE)dl + c^~f3{E)^2TTn. (1) 
Jl Jr 

Here, the first and second terms represent the phase gain 
along the trajectory, (f> is the phase difference between 
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FIG. 1. Grain boundary junction between two d-wave su- 
perconductors. Here a and 6 denote the crystalline axes, and 
5x is the mismatch angle. The y-sxis is chosen along the 
grain boundary. The diagrams in the middle indicate the 
directions connecting the lobes of the order parameter with 
the same sign (shaded) and opposite sign, on both sides of 
the grain boundary. Supercurrent in the shaded directions is 
carried by regular Andreev levels, otherwise by 7r-levels (see 
text). For given Sx, the properties of the junction depend 
on the grain boundary orientation vs. axes a, h. The special 
cases are asymmetric (xi = 0) junction, (b), and symmetric 
junction, xi = —Xr (c). The positive and negative lobes of 
the order parameter are chosen along a and b respectively (d). 



the order parameters in the points L, R (including the 
intrinsic phase difference), and the additional phase shift 
(i{E) = IT ■ 0{1) depends on the shape of the scattering 
potential A(r) near the reflection points (/3(0) = tt). The 
minus sign before the second term reflects the quasipar- 
ticle branch change. Evidently, the right-moving elec- 
tron and left-moving hole both carry electric current in 
the same direction. Therefore, unlike conventional stand- 
ing wave solutions to the Schrodinger equation, Andreev 
bound states (Andreev levels) (Eq. (|l])) carry electric cur- 
rent. This provides a mechanism for the Josephson effect 
in such structures as microbridges or in SNS contacts 
|^t|-|13|, which can be extended to the case of a general 
Josephson junction (see |Q). The current through a 
junction of unit width is expressed as pq| 



ek^ ^ r'"^ 
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Here, the a;-axis is chosen normal to the grain boundary 
(see Fig. |l|), np{E) is the Fermi distribution function, and 
En{4>,0) is the energy of the n-th Andreev level for an 
electron with incidence angle 9. The direction-dependent 
intrinsic phase of A in d-wave superconductors leads to 
qualitatively new features of surface/interface states in 
d-wave superconductors (for current reviews see [ |l6| , p^ ). 
They stem from the coexistence of two types of Andreev 
levels, which are formed by Andreev reflections from the 
order parameter with the same (regular levels) or differ- 
ent (tt- levels) intrinsic phases (see Fig. |l|). In case of a 
d-wave superconductor boundary with an insulator, the 



latter include so called midgap states (MGS) , or zero en- 
ergy states (ZES) exactly at the Fermi level (since 
for (/) = TT the condition (|^) is always satisfied by i? = 0) . 
The two sets of Andreev levels carry Josephson currents 
in opposite directions. In clean SND or DND junctions, 
as well as in short SD or DD point contacts, their con- 
tributions behave otherwise similarly. As a result, the 
current-phase dependence can become 7r-periodic [|l9|-pi| . 
What is more interesting, the equilibrium phase differ- 
ence across such a junction is neither 0, nor tt, but is 
degenerate, ±0o, which reflects the T-breaking ground 
state of the system 0. The value of (/)o at T = de- 
pends only on the orientation of the d-wave order pa- 
rameter with respect to the boundary pT[ |. The occur- 
rence of spontaneous currents in this picture is natural 0| 
(Fig. In equilibrium, the Josephson current between 
the left and right superconductors (in the x-direction) is 
zero, which is possible only if the contributions of regular 
and TT-levels cancel each other. But the y-components of 
these contributions add up, yielding a spontaneous cur- 
rent along the boundary. Evidently, there are two equi- 
librium states (with a spontaneous current flowing up or 
down). 




FIG. 2. Andreev levels and spontaneous currents in a DND 
model. The total superconducting current across the bound- 
ary, carried by normal levels (a) and 7r-levels (b), is zero in 
equilibrium, while the spontaneous current js parallel to the 
boundary is finite. The degeneracy of the ground state is 
illustrated by the freedom of choice of its direction. 

This picture becomes more complicated in the pres- 
ence of finite boundary transparency D, which breaks the 
symmetry between zero and 7r-states. The latter are not 
split by specular tunneling, and as a result, the ampli- 
tude of supercurrent carried by them scales as y/D due 
to resonant transmission p5| ] (instead of D for regular 
levels) . This led to the prediction of temperature-driven 
transition to 7r-junction [ p2| in a DD junction. Recent 
measurements of 7(0) in YBCO grain boundary junctions 
p3| generally confirmed this prediction, but showed that 
the transition does not necessarily occur between equi- 
librium values (l)Q = and (po = it. 

Notice that in all of the above arguments we did not as- 
sume any subdominant order parameter. In other words, 
the T-breaking and generation of the spontaneous cur- 
rent happens without any subdominant pairing potential. 
In general, a subdominant order parameter will occur 
at the boundary, if there exist interactions in the cor- 
responding channel. In DD-junctions, the source of the 
subdominant order parameter is the proximity to a differ- 
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ently oriented superconductor and therefore it will exist 
at all temperatures even above the subdominant critical 
temperature Tc2- It is interesting that the presence of a 
subdominant interaction channel will actually suppress 
the spontaneous current pj] . Wc shall discuss this coun- 



F(vi?,v^) [see Eq. (A4)]. In a pure dx2_y2 case we as- 
sume the interaction potential to have the form 



terintuitive behavior in section III 



The role of the subdominant order parameter in the 
spontaneous current generation at surfaces is completely 
different. Near a boundary [e.g. the (llO)-surface] the 
dominant d^2_y2 order parameter is suppressed |p5|| and 
therefore a subdominant order parameter {d^y or s) can 
appear. It will be formed below the smaller critical tem- 
perature Tc2 if there exists nonzero pairing interaction in 
the corresponding channel |^. The combination of the 
two order parameters with complex coefficients breaks 
the T-symmetry [Q| and can lead to spontaneous sur- 
face currents and magnetic fluxes. Usually d.j.2_y2 ± is 
or dx^^y^ ± idxy combinations are predicted. Recent ob- 
servations of zero bias peak splitting in surface tunneling 
experiments [^ and spontaneous fractional flux (0.1-0.2 
<I>o) near the "green phase" inclusions in YBCO films ||2^ ] 
agree with this picture. 

The simple description based on Andreev levels pre- 
sented earlier is best suited for quasi-lD problems in SNS 
(or DND) structures. It is possible to generalize it to 
deal wit h at l east some of the above mentioned compli- 
by introducing "Andreev tubes" of width 



cations |29 



^ Xf, following the quasiclassical trajectories. Neverthe- 
less, in order to obtain quantitative results, it is better to 
use an approach based directly on the method of quasi- 
classical superconducting Green's functions (Eilenberger 
equations [3l|). In this paper we apply the formalism 
to the case of a planar d-wave superconductor in contact 
with another superconductor or vacuum, for an arbitrary 
transparency and roughness of the boundary. 

The paper is organized as follows. In Section ||, we 
consider the generation of spontaneous currents by the 
proximity effect at a uniform SD or DD interface. We 
obtain the equilibrium phase, the spontaneous current 
distribution, and the superscreening effect (in the latter 
case) in case of ideal, rough, and partly reflective surface. 
In Section III , the spontaneous current generated by the 



subdominant order parameter at a boundary is consid- 
ered. The interplay of the two mechanisms (proximity 
effect and subdominant pairing) is discussed for DD and 
SD junctions. The technical details of the formalism are 
given in the appendices. 

II. SPONTANEOUS CURRENT GENERATED BY 
PROXIMITY EFFECT 

A. Ideal junctions 

Let us consider a planar d-wave superconductor with 
a straight grain boundary along the y-axis in its ab- 
plane (cf. Fig. |l|). The order parameter A(vi?,r) is 
self-consistently determined by the interaction potential 



VivF,Vp)^VdCos2{e-x)cos2{6' -x), 



(2) 



where the angles 9, 6' give the direction of v^?, in the 
a6-plane, and x is the angle between the crystallographic 
a-direction and the z-axis. The dimensionless BCS con- 
stant of interaction is — VdN{Q)/2. We can also con- 
sider a boundary between a d-wave superconductor and 
an s-wave film in which case, for the s-wave supercon- 
ductor, we have 



(3) 



This problem is essentially 1-dimensional, with 
A(vi?,x) — > Ar,i(vi?) as X — > ±oo, where the subscripts 
I and r represent left and right of the boundary respec- 
tively. The method we choose to solve this problem is 
the standard quasiclassical method which is described in 
detail in appendix |^. The self-consistent solution can be 
obtained only numerically (see appendix |c|) . However, 
we start from the analytical solution in the simplest ap- 
proximation, when we assume that the order parameter 
takes its bulk values at all x, there is no subdominant 
order parameter, and the grain boundary is ideal (trans- 
parent and specular, see appendix P). 

We denote by jj [= jx{x=Q)\ the Josephson current 
flowing from the left to the right superconductor, and by 
is [= jy{x—Qi)\ the surface current flowing along the in- 
terface at the boundary. All the current distributions are 
expressed in units of [defined in (A6)] which is of the 



order of the bulk critical current density. These currents 
are expressed by Eqs. (B8) and (B9) which are valid 
(within the applicability of the model) for arbitrary sym- 
metry of the order parame ters A;.,.. For a DD interface, 
the functions Ai^r(vF) in (p8|) and (|B9[) are 



A;,, = Ao(r) COs2(0-Xi,r), 



(4) 



where An(r) is the maximum gap, as introduced in ap- 
pendix |D|. Note that in this section the pairing potential 
is assumed to be nonzero only in a single orbital chan- 
nel on either side of the boundary. Therefore, while the 
anomalous Green's functions (/, /^) which support differ- 
ent orbital symmetries are induced across the boundary, 
they do not translate into a subdominant order parame- 
ter A'. 

The results of the calculations for a DD junction are 
displayed in Fig. |^ for different mismatch angles between 
the crystalline axes across the grain boundary and at 
temperature T = O.lTc (assuming the same transition 
temperature on both sides). In all these figures, the left 
superconductor is assumed to be aligned with the bound- 
ary while the orientation of the right superconductor 
varies. The Josephson current-phase relation (Fig. ||a) 
demonstrates a continuous transition from a 7r-periodic 
(sawtooth-like) line-shape at 5x = 45° to a 27r-periodic 
one for small 5x, as expected in the case of a clean DND 
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FIG. 3. Josephson current (a) and spontaneous current (b) 
versus the phase difference in a clean DD grain boundary junc- 
tion calculated in non-self-consistent approximation . Current 
densities are in units of jc [cf. Eq. (A.6)] and T = O.lTc- The 
mismatch angles are Xi = and Xr = 45° (1), 40° (2), 34° 
(3), 22.5° (4). 



junction pT| . The phase dependence of the surface cur- 
rent (Fig. Mi) is also in qualitative agreement with earlier 
results for SND and DND junctions 0. 

It is important to make a cautionary remark here. 
In order to have the Josephson effect, there must be a 
weak link between the two superconductors. In other 
words, the superconducting phase should change over a 
short distance at the boundary (otherwise one should 
speak about phase gradient and not phase difference). 
The weak link in a clean DND junction is provided by 
the normal layer, while in an ideal DD junction, it is 
formed due to the suppression of the order parameter 
near the boundary (at finite 5x) , which follows from self- 
consistent treatment. At Sx = 0, such a "junction" (of 
infinite width) is simply a bulk superconductor and not 
a weak link, therefore Josephson physics does not ap- 
ply. Nevertheless, the non-self-consistent approximation 
is applicable, in the limit Sx — > 0, to the case of a narrow 
contact (microbridge) between the two sides, in a way 
similar to conventional superconductors p3]. 
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FIG. 4. Equilibrium phase difference in a clean DD grain 
boundary junction as a function of 5x = Xr — Xi (keep- 
ing xi =0), at different values of t = T/Tc- The circles 
and triangles correspond to non-self-consistent (NSC) and 
self-consistent (SC) calculations respectively. For nonzero (po, 
the ground state is twice degenerate (0 — ±00 ). 



The equilibrium phase difference across the junction, 
at which jj{4>) = and dj j {(j)) / d(j) > 0, takes any value 
between and 7r/2, and is degenerate, (f> = ±4)0, unless 
(j)o = 0. The T-symmetry is therefore broken. The re- 
gion of T-breaking states (as a function of temperature 
and mismatch angle) is shown in Fig. ^. In that fig- 
ure we also present the self-consistent numerical result 
for comparison. The method adopted for our numerical 
calculations is described in detail in appendix |^. Only 
in the asymmetric Sx = 45° junction does the degener- 
acy (at (j) = i'''/2) survive at all temperatures, due to 
its special symmetry which leads to complete suppression 
of all odd harmonics of /(0); generally, (/)o — > at some 
temperature that depends on the orientation. The equi- 
librium value of the spontaneous current is nonzero in a 
certain region of angles and temperatures (Fig. |5|), which 
is largest in the case of asymmetric Sx = 45° junction. 
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FIG. 5. Spontaneous current in the junction of Fig. ^. 

We can also calculate the function jy (r ) at all x us- 
ing our analytical formalism [cf. Eq. ( BIO )]. This func- 
tion is plottedrin Fig. In this figure we also plot the 
same graph calculated using the self-consistent numeri- 
cal method described in appendix with and without 
the subdominant order parameter (subdominant pairing 
is discussed in appendix ^ . The curves are qualitatively 
similar although there is a small quantitative difference. 
The same is true for the dominant order parameter, but 
not for the subdominant order parameter which appears 
only when the interaction in the corresponding channel 
exists (cf. Fig. A in this figure and all subsequent 
figures is normalized to Tc). Notice that the subdomi- 
nant order parameter exists near the boundary despite 
T > Tc2- The reason is that the appearance of the 
subdominant order parameter is merely the result of the 
proximity effect, i.e., the extension of the order param- 
eter from one region to the other, which happens at all 
temperatures. 

In all of the graphs in Fig. ^ we see a remarkable fea- 
ture: the current along the boundary is sharply peaked in 
a layer of order around it, but is accompanied by coun- 
terflows spread over about lO^o on either side. Within 
the numerical accuracy, the total current in y-direction 
is zero, independently on either side of the junction. 
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FIG. 6. Equilibrium spontaneous current distribution (a) 
and order parameter (b) for a (0°-45°) grain boundary junc- 
tion at t = 0.1. Solid and dashed lines represent the re- 
sults of non-self-consistent and self-consistent calculations 
respectively. Long-dashed and dot-dashed lines represent 
the case with interac 
(tc2 = Tcz/T; = 0.05). 



Since in high- Tc superconductors lO^o <C Al,Aj (Lon- 
don and Josephson penetration depth res pec tively), the 
phenomenon can be called super screening [p4| . Note that 
we so far did not take into account the Meissner screening 
- and as it turns out don't need to; the magnetic field of 
spontaneous currents is cancelled on a scale less than Al • 
This effect may be responsible for difficulties with ob- 
serving spontaneous currents at surfaces and interfaces 
of d-wave superconductors [Q, although as we will see, 
the effect is suppressed by roughness or reflectivity of the 
boundary. 
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FIG. 7. (a) Spatial distribution of spontaneous current 
density in SD and asymmetric DD clean junctions at 
t = tc2 = 0.05 and tcs = 0.1. In both cases xi ~ 45° and in 
the DD case Xr ~ 0. (b) Spontaneous current near DD junc- 
tions with different misorientation angles Xi = and Xr = Sx 
at i = 0.1. Calculations in this figure and in all subsequent 
figures are done self-consistently. 

In the case 6x = 45°, the superscreening can be ob- 
tained analytically from Eq. (BIO); the integral 



uj>Q 



Al Ar sin 6 sign (cos 6) sin \vp cos 8 \ 



P- + A; Ar cos ( 



is exactly zero after angular averaging. However, self- 
consistent numerical calculations show (within the nu- 
merical accuracy) the same behavior at all orientations 
(Fig. 0b). To better understand the situation, let us 
recall that in a system with local magnetic moment 
density m(r) the "molecular currents" flow with den- 
sity j(r) =cV X m(r) (Fig. ||). In a superconductor 
with order parameter d,j.i_yi -I- e^'^^d^y, the local or- 
bital/magnetic moment density is 

/'27r in 

l(r) cx m(r) cx z / — (Ai(a;) cos 26* + A2(a;)e"''^« sin 26l) 

x--^(Ai (x) cos 29 + A2 (x)e''l"> sin 20) 
I oO 

= 2Ai(x)A2(a;)z sin^o- 
The contribution to the spontaneous current is thus 
j(r) oc V X l(r) || y. Note that the same expression is 
obtained from the Ginzburg-Landau equations p4[ : 
j oc V X z Im[(ii (r) (jj (r)]- The total current in 
the y-direction due to this mechanism is Itot oc 
dS • V X m cx ds ■ m, where fi is a cross-section 
in the xz-plane (assuming the superconducting film has 
finite thickness h), and ds is the linear element of its 
boundary dfl. The latter integral is obviously zero, be- 
cause either m || z _Lds (when z = zLh/2), or m = 
(at X = ±00 , where the contour dfl is closed). More- 
over, the "molecular currents" picture, shown in Fig. |^, 
reproduces to the distinct peak-and-counterflow current 
distribution. 



\ 

qOQ 



\ 

QQo 



dx jy{x) oc 



(5) 



FIG. 8. "Molecular currents": the supercurrent distribu- 
tion in a DD junction as the superposition of local currents 
produced by the angular momenta of Cooper pairs in a chiral 
d + e^^'^«d' state. 

A finite magnetic moment in the bulk is essential for 
the effect. The corresponding chiral state is created by 
the proximity effect, when the two d-wave order param- 
eters mixed near the boundary are spatially rotated by 
a nontrivial angle 5x and have a nontrivial phase shift 
±(/)o 7^ 0, IT. Although the argument given here is based 
on having a mixed order parameter near the boundary re- 
gion, in fact the presence of the second order parameter is 
not necessary. The reason is that the subdominant sym- 
metry already exists in the Green's functions g^j and f^i, 
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even if the interaction potential does not support pairing 
in that channel. Since the current density is related to 
these functions rather than to the pairing potential [cf. 
Eq. (A5)], the molecular current is basically formed by 
inherently mixed symmetries of these functions. 

The above arguments are certainly not true when the 
order parameter symmetry is c? + is, as in SD junctions 
(Fig. 0a). In this case the molecular current is identically 
zero and as a result no countercurrent exists. (Of course, 
since the Meissner currents must be taken into account 
in this case, the results are valid only if the system size 
is much less than the London penetration depth.) 



B. Reflective junctions 

We also study the properties of the spontaneous cur- 
rent in the presence of a non-ideal boundary, i.e., with 
nonunity transparency and/or surface roughness. The 
details of our method of numerical calculation are given 
in appendix The transparency and roughness of 
the boundary are parameterized by < Dq < 1 and 
< P < oo respectively (see appendix ^ for definitions). 
In particular, in the case of a clean ideal junction, Dq = 1 
and p = 0. 
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FIG. 9. Comparison between the size of the spontaneous 
current in asymmetric and symmetric Sx = 45° DD junctions. 
Both calculations are done at t = 0.5. 

An interesting case is presented by a symmetric 45°- 
junction (x; = -22.5° and Xr = 22.5°). Although the 
ground state is degenerate in this case pS], the sponta- 
neous current is practically absent (Fig. This is easy 
to see from geometric considerations (Fig. 0c). Consider 
the total supercurrent across the boundary as a sum of 
contributions from quasiclassical trajectories. There are 
two groups, analogous to Andreev regular and 7r-levels 
discussed in the introduction; the corresponding direc- 
tions are within the shaded and white sectors in Fig. |I]c 
respectively. In equilibrium, their contributions to the 
current normal to the boundary must cancel each other. 
What happens to the tangential component of the equi- 
librium (spontaneous) current, depends on the orienta- 
tions of the order parameters with respect to the bound- 
ary and each other, and the equilibrium phase difference 
(also a function of temperature). It is obvious from the 
picture that for any given mismatch angle (5x, in a sym- 
metric junction, the tangential components of the cur- 
rent cancel for regular and 7r-directions separately. This 
is why symmetric junctions are very "quiet" (using the 
term of 139]): they can violate T-symmetry without pro- 
ducing local magnetic fields, which could couple to some 
external degrees of freedom 40 1. This can be both an 
advantage and a disadvantage from the point of view of 
using such a system as a solid-state qubit |39,pq]. 
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FIG. 10. Order parameter (a) and spontaneous current dis- 
tribution (b) in (0°- 45°) asymmetric grain boundary junc- 
tions with different transparencies Do at t = 0.1. 
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FIG. 11. Influence of finite barrier transparency on sponta- 
neous current amplitude in an asymmetric (0°-45°) DD junc- 
tion without (a) and with (b) subdominant order parameter. 
Here "L"= left and "R"= right of the boundary and the cal- 
culations are done nt t = 0.05 and tc2 ~ 0.1 (for (b)). In 
the case of (b), one can see that as Dq —> 0, the spontaneous 
currents on the right hand side of the boundary survive. 



In Fig. |10| a we see the self-consistent pairing potential 
near the interface for several values of Dq. Notice that 
for the case of Dq — (i.e., totally reflecting surface), 
the order parameter is constant on the left while it van- 
ishes at the junction on the right side of the boundary. 
The latter is because the quasiparticles on the right see 
opposite sign of the order parameter after reflection from 
the boundary As Dq — > 1, the order parameter 

becomes continuous at the boundary. The spontaneous 
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current (Fig. |lO|b) also evolves from zero to its maxi- 
mum value as Dq changes from to 1. Notice that at 
all values of Dq except for and 1, the magnitudes of 
the spontaneous current on the left and right hand sides 
of the boundary are different and exact supcrscreening 
happens only in the ideal boundary case {Dq = !)• Fig. 
|TT| a demonstrates the continuous evolution of the spon- 
taneous current from zero to its maximum as a function 
of Dq. This picture is influenced by the presence of a 
subdominant interaction in Fig. Especially, on the 

right hand side of the junction (with Xi- — 45° orienta- 
tion), the spontaneous current does not vanish even at 
Do = 0. We will come back to this point in the next 
section. 

C. Rough junctions 

The effect of surface roughness is very different from 
that of finite transparency. First notice from Fig. |l^ that 
surface roughness has little effect on the order parameter 
distribution (Fig. [T^a), while it significantly changes the 



spontaneous current (Fig. |12p) . This may seem strange 



until we recall that supercurrents (including spontaneous 
currents) are generated not by the order parameter A, 
but by the Green's funct ions , /^^ and g^^ [related by the 
normalization condition (A3)]. The latter is directly af- 
fected by the boundary roughness, while the former isn't. 
Again, we notice that exact superscrccning is absent at 
finite roughness. 




-4 -2 2 4 6 

FIG. 12. Order parameter (a) and spontaneous current dis- 
tribution (b) in (0°-45°) asymmetric grain boundary junction 
witli different values of the surface roughness p at f = 0.1. 

Another interesting effect is presented in Fig. |lja: as 
the surface becomes rougher and rougher (p oo), the 
spontaneous current on the left side (with xi = 0) van- 
ishes while on the right side {xr = 45°) it saturates to 
a finite value. This non-trivial behavior is directly re- 
lated to the anomalous proximity effect between a d-wave 
superconductor and a disordered region studied by Gol- 
ubov and Kupriyanov js^ . When p is large, the contribu- 
tions from different quasiclassical trajectories are mixed. 
Therefore the angle-dependent components of /, which 
"leak" through the boundary to the other side, are sup- 
pressed, but in the limit of infinite roughness the s-wave 



contribution still survives. The latter is generated by 
scattering from the rough boundary even in case of only 
dominant d-wave pairing on either side, and can be es- 
timated as fs,eB = J^l^/^{de/Tr)cos9f{e,0^) [the upper 
(lower) sign corresponds to leakage from the right (left) 
to the left (right) side of the junction]. In an asymmet- 
ric (0°-45°) DD junction we thus expect exactly zero in- 
duced / on the left due to symmetry, while on the right 
a finite s-wave component should appear with amplitude 
2/{3tt) w 0.2 relative to the d-wave anomalous Green's 
function across the boundary. The spontaneous currents 
in this case can flow only on the right of the boundary, 
where the chiral combination d ± is is thus formed, in 
agreement with our numerical results. This might also 
explain why supcrscreening is absent in case of rough 
surfaces: d + is symmetry does not supprot molecular 
moments. Non-ideal asymmetric grain boundaries be- 
tween d-wave superconductors can thus be better candi- 
dates for the search of spontaneous currents, since the 
supcrscreening is suppressed, while the supercurrent am- 
plitudes are still detectable. 
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FIG. 13. Influence of surface/interface roughness on the 
magnitude of spontaneous current in a junction (a) and near a 
surface (b). Note that in the latter case spontaneous currents 
are suppressed at much smaller values of p. All calculations 
are done at t = 0.1 



III. SPONTANEOUS CURRENT GENERATED 
BY SUBDOMINANT ORDER PARAMETER 

We already saw the effect of the presence of a sub- 
dominant order parameter on the spontaneous current in 
DD grain boundary junctions (Fig. ||). In particular in 
Fig. ^ we notice that the subdominant component ex- 
ists even at temperatures above the subdominant critical 
temperature, Tc2- The subdominant order parameter in 
Fig. 1^ is therefore completely induced by the proxim- 
ity to a differently oriented superconductor on the other 
side of the junction. On the other hand, in Fig. pl] b we 
realized that the spontaneous current exists even when 
Dq 0, i.e., when the proximity effect is completely ab- 
sent. This happens only when T < Tc2 and therefore the 
mechanism is completely different from the one described 
above. 
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In general, the bulk of a superconductor can support 
only one symmetry of the order parameter even if the 
interaction potential contains finite interaction in several 
different channels (unless the coupling constants of the 
channels are very close to each other ^^). This is be- 
cause the dominant order parameter introduces a cutoff, 
in the BCS gap equation, that removes the logarithmic 
divergence (at small w), responsible for the exponential 
dependence of the gap on the interaction potential (see 
appendix ^ . As we saw in Fig. |l^a, at a reflective surface 
of a 45° oriented d-wave superconductor, the dominant 
order parameter vanishes. Suppression of the dominant 
order parameter gives the chance to the next subdomi- 
nant one to appear if the temperature is below its cor- 
responding critical temperature. Thus, the appearance 
of subdominant order involves a spontaneous symmetry 
breaking and therefore a second order phase transition. 
This is demonstrated in Fig. |l^ (we increase Tc2 in- 
stead of T in that figure, but the behavior is the same). 
The subdominant order parameter always appears with a 
phase 7r/2 with respect to the main one. With this phase 
difference, the mixed order parameter (d^2_y2 ±idxy or 
dx2-y2 ± is) is fully gapped (with no nodes), which is 
energetically favorable because of the extra condensation 
energy gained. 

It is important to notice that the above phenomenon is 
not likely to happen at the interface between supercon- 
ductors because the dominant order parameter is never 
completely suppressed (except when Z) = 0, cf. Fig. 10 1). 
Even at rough interfaces, the dominant order parameter 
at the boundary is only slightly less than half of its bulk 
value (Fig. [T2|a). The subdominant order parameter in 
this case is predominantly induced by the proximity ef- 
fect and therefore follows the symmetry of the neighbor- 
ing superconductor. This may explain why the mixed 
symmetry state was not observed in a recent experiment 



conductor are presented in Fig. |14|a. Qualitatively, they 
are the same as at a DD or DS boundary, with the same 
superscreening behavior in the case of d -I- id' order pa- 
rameter symmetry. Despite the similarities in the current 
distribution, there are fundamental differences between 
the two cases. First of all, unlike at the interface, the 
appearance of spontaneous current near the surface is a 
very fragile phenomenon. A small deviation from the 45° 
angle will suppress the effect significantly. It is also very 
sensitive to surface roughness as compared to the inter- 
face current; a small p is enough to remove the spon- 
taneous current (Fig. [l^b). Moreover, unless the effects 
of the magnetic field generated by spontaneous currents 
itself are taken into account, the presence of a subdomi- 
nant gap is necessary for the surface spontaneous current, 
but not for the interface one. (The spontaneous symme- 
try breaking due to Doppler shifts of midgap states ^] 
occurs at temperatures below 1/6A(^o/Al) ji^-^^, that 
is approximately IK for high-T^ compounds.) 



(a) DD Junction 



(b) SD Junction 




FIG. 15. Suppression of spontaneous currents in the pres- 
ence of subdominant order parameter as a function of sub- 
dominant transition temperature. In both cases we have 
taken the same Tc on both side of the junction. 



0.05 




FIG. 14. (a) Spontaneous current profile near the 
(llO)-surface of a d-wave superconductor with d^y- and 
s-wave subdominant order parameter. Here t — 0.05 and 
tc2 ~ 0.1 for both the s and dxy subdominant gaps, (b) 
Spontaneous current at the surface as a function of the sub- 
dominant transition temperature at t = 0.1. 

Spontaneous currents at a surface of a d-wave super- 



It is interesting to note that the presence of interaction 
in a subdominant channel actually works against the ap- 
pearance of spontaneous current at an interface . This 
behavior is displayed for both DD and SD junctions in 
Fig. |l5|. The reason becomes clear from considering the 
Andreev levels in the junction (which can be modeled on 
this occasion by a DND junction. Fig. ^6|). The induced 
subdominant order parameter will be aligned with the 
dominant one across the boundary. Therefore in addition 
to the "dominant-dominant" set of Andreev levels there 
will appear a "subdominant-subdominant" set, which ob- 
viously carries supercurrents in the opposite direction. Of 
course, a subdominant pairing potential should also sup- 
pres s the Josephson current, which was indeed noted in 

As we have mentioned above, the subdominant order 
parameter near interfaces is induced across the boundary 
from the other side, where it is dominant, and can be sup- 
pressed only simultaneously with the latter. Therefore, it 
seems impossible to directly observe the influence of the 
subdominant order on the spontaneous currents by e.g. 
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FIG. 16. The physical mechanism of spontaneous current 
suppression by subdominant order parameter in DND model 
(see text). 
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APPENDIX A: QUASICLASSICAL 
EILENBERGER EQUATIONS 



raising the temperature above Tc2 to suppress it. The 
same, of course, holds true for suppression by magnetic 
field. 



IV. CONCLUSIONS 

We have investigated the occurrence of spontaneous 
currents at surfaces and interfaces of d-wave supercon- 
ductors, using the self-consistent quasiclassical Eilen- 
berger equations. A T-breaking ground state is the 
only necessary condition for their appearance, with no 
quantization condition for the generated magnetic fluxes. 
Therefore, the effect is sensitive to the properties of the 
system and allows the existence of "quiet" T-breaking 
states, i.e., states with fluxes much smaller than the flux 
quantum. 

We have shown that the spontaneous current at 
the (llO)-surface of a d-wave superconductor is formed 
through a second order phase transition to a mixed sym- 
metry state (d -|- id! or d -\- is), and is very sensitive 
to the interaction in the subdominant channel, temper- 
ature, and surface roughness and reflectivity. On the 
other hand, at an interface between two differently ori- 
ented d-wave superconductors, or between a d-wave and 
an s-wave superconductor, the spontaneous current is 
generated as a result of the proximity effect. It is gener- 
ally robust, with less sensitivity to the above effects. In 
particular, at very rough surfaces, the spontaneous cur- 
rent survives on one side of the junction, in agreement 
with the recently proposed anomalous proximity effect 
at rough d-wave surfaces. We also show that interaction 
in any subdominant channel suppresses the spontaneous 
current at the interface, while its existence is necessary 
for the T-breaking at the surfaces. 

The contribution to the spontaneous current from the 
local orbital moment of the condensate leads to char- 
acteristic "superscreening" of the equilibrium current in 
DD junctions and at surfaces with dxy subdominant or- 
der. Exact superscreening is sensitive to boundary im- 
perfections and disappears at flnite roughness and/or 
nonunity transparency of the junction. 



To describe the coherent current states in a supercon- 
ducting ballistic microstructure we use the Eilenberger 
equations for the ^-integrated Green's functions ]3l| ]: 

Vi=- ■ Q^G^ivF, r) + [ljts + A(vF, r), G^{vf, r)] = 0, 

(Al) 



where 



A = 



A 
At 



9ui fij 



(A2) 



A is the superconducting order parameter and G^ji'vp, r) 
is the matrix Green's function, which depends on the 
electron velocity on the Fermi surface vp, the coordinate 
r and the Matsubara frequency ui = (2n -I- l)7rr, with 
n being an integer number and T the temperature. We 
also need to satisfy the normalization condition 



(A3) 



In general, A depends on the direction of vj;^ and is de- 
termined by the self-consistency equation 

A(vi., r) = 27rA^(0)T ^ (F(vj^, v^)/^(v^, r))^^ (A4) 

a)>0 

where V{vp, v'p) is the interaction potential. Solution of 
matrix equation (Al) together with (A_4) determines the 
current density j(r) in the system 



j(r) = -AmeN{0)Tj2 i^FgU^F^r))^ 



(A5) 



w>0 



In two dimensions, A^(0) — m/27r is the 2D density of 
states and (...) = / {dO /2tt)... is the averaging over direc- 



tions of 2D vector v^?. Throughout this article we write 
current densities in units of 

= ATT\e\vFN{Q)T, (A6) 

which is of the order of the bulk critical current density. 
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Supposing A(— vp) = A( v^), which is always the case 
for superconduc tor s with sing let p airing, from the equa- 
tion of motion ( Al) and Eq. (A3) we have the following 
symmetry relations: 



/*(w,-VF 

5*(w,"VF 

g{-uj,-VF 
At 



.9(w,Vi^), 
/(^,vf), 

A*. 



APPENDIX B: NON-SELF-CONSISTENT 
ANALYTICAL SOLUTION 

We assume a boundary, at a; = 0, between two super- 
conductors. To find an an alyt ical solution we neglect the 
self-consistency equation (A4), and write 



A(v^, .) = ( " ; n ' (Bl) 

^ ^ [ Ar(vi?) exp(— 10/2), a; > ^ ' 

where I (r) stands for left (right) of the boundary. Ai^r 
can in general have d-wave, s-wave, or any other symme- 
tries. The Eilenberger equations for the Green's function 
Guj are linear and therefore can be easily solved sepa- 
rately for positive and negative x to satisfy 



A; r <^ 

lim = — ^, lim = -— , 



(B2) 



where fli^r — Y^^ + I^i,r| 
solutions. For x < 



This yields the following 



for a; > 



(7/f7; -cj)e2^"'/l"-ICi 



„ +^_(-^17,-c.)e2-"'/l'^^ICi 



(B3) 



^ -I- ^^{-V^^r - c.)e-2-"'Vk.lc2 

A^e*'^/^ 



J<h/2 



(B4) 



Here 77 = sign(wa;). Imposing the continuity condition for 
Guj at the boundary we find 



C2 



Al uj(Ai — Ar cos (/)) — irjAr^i sin q 

fll iliilr + Uj'^ + AiAr cos (j) 

Ar uj{Ar — Al COS0) — irjAiflr sine 

fir ftlflr + OJ^ + AlAr COS 6 



(B5) 



Substituting into (B3) or (B4) at a: = we find the 
Green's functions / and g to be 

/(O) = 00 , ^ .2 , A A ^ ' (B6) 



ff(0) = 



QiQr + io^ + A/ Ar cos ( 
ll}{Qi + fir) ~ i??A;Ar sin0 

fllflr + Oj'^ + A; Ar cos (j) 



(B7) 



Using the expressions (A4) and (B7), we obtain the 
current densities jx{0) = j.i and jy{0) = js 



JJ 



u)>0 



UJ>0 



AiAr Icos ( 



fllflr +0^2 + A; Ar COS ( 



A; Aj. sin 6 sign(cos 0) 

fllflr + ijj'^ + A; Ar COS ( 



sm( 



(B8) 



sin 0, (B9) 



whe re the current densities are in units of jc [defined in 
(|A(]|)] and t = T/Tc- At all other points we can find j 
from 



j(^)-tE 



vp 7/A/Ar sin <f> e 



-2\x\ilj^ / \v F COS 6\ 



fllflr -1-07^-1- A; Ar COS ( 



(BIO) 



where vp is a unit vector in the direction oi wp. 

APPENDIX C: SELF-CONSISTENT 
NUMERICAL SOLUTION 

The general approach we use is based on transforma- 
tion of the set of coupled (formally) linear Eilenberger 
equations (Al) for the functions f,f\g to two nonlin- 
ear, but separate, equations which are numerically sta- 
ble (Schopohl-Maki transformation |^^). To this end we 
express the components of G^^ matrix as 



1-ab , 2a . 2b 

9^1-. — r. 1 = 1-. — r> / = 



\ + ah 



l + ah 



1 + ab 



(CI) 



Now a and h satisfy two independent nonlinear equations 

vj;^ • Va = A — A* — 2uja 
-MP ■Vb = A* - Ab^ - 2ujb. (C2) 

Notice from these equations that a{—\-p) = b*{\-p) and 
b{—vp) — a*{vp). We use the solutions for a homoge- 
neous system, 



A 



(C3) 



where fl — ^/uJ^+\A\^ , as asymptotic conditions at 
X ±00. For positive Vx, the first (second) of Eqs. 
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(|C2|) is stable if we choose the boundary condition at 
— oo (+CX)). Using the appropriate asymptotic condition 
one can find a (b) at aU other points on the quasiclassi- 
cal trajectory along the vector v^, integrating along the 
trajectory. Self-consisten cy i s achieved iteratively. The 
step- like approximation (Bl) is used to find a, 5 in the 
first approximation, which are in turn substituted in Eq. 
(A4) to find the next iteration for A. These steps are 
repeated until A does not change within numerical accu- 
racy. 

It is convenient to take the order parameter constant 
between discrete points on the trajectory (the angle 9 
gives the direction of v^? as usual), separated by a dis- 
tance h. Then we find 



A, - c^A* 



2aiU} 



M+i 



bi-i = bl 



QiA* + UJ - 

A* - 



- ili coth{rtih/ cost 
6? A, - 2biLU 



biAi + UJ + coth(rii/i/ cos 9) 



(C4) 



(We have explicitly taken into account that this proce- 
dure is stable in opposite directions for a,b). Having 
obtained a and b, we can find / and g and therefore the 
orde r parameter A and current density j using ( A4 ) and 



1. Effect of surface reflectivity 

Our numerical method can also be applied when the 
transparency of the junction is arbitrary, < D < 1. 
Since part of the quasiparticles get reflected from the 
boundary, the quasiparticle trajectories of the reflected 
quasiparticles and the transmitted ones from the other 
side will mix. Then even in the no n-se lf-consistent ap- 
proximation for the order parameter (Bl) we cannot sim- 
ply impose continuity of the Green's functions along a 
trajectory as we did before for the ideal transparency 
case. Instead, one should use Zaitsev's boundary condi- 
tions which for > is given by 



D 



D 




- 12 



(C5) 



(C6) 



where 



(C7) 



with being the reflection of vj? with respect to the 
boundary. Similar relations also hold for s ' and d ' . For 
Vx < one has to replace r <-^- Z in ( p6| ). In general, 
D can be momentum dependent. In our calculations we 
consider [nol 



Dn cos^ 9 



1 - Do sin^ 9 
Let us introduce four Green's functions 



Gi = G^(0_,( 
G2 = G(j(0+, ( 
We can now write 



G3 = G(^(0_,7r - 

G4 = Gcj(0+, TT — I 



(C8) 



(C9) 



s' = Gi + G3, 



Gi — G3 
G4, d ^ = G2 ^ G4 



(CIO) 



Rewriting these Green's functions in terms of ,bi, {i^ — 
1,...,4), we can calculate 01,62,63, and 04 by integrat- 
ing from the corresponding infi nity in the direction of 
convergence [e.g. using Eq. (C4)]. At the boundary, the 
remaining functions 61 , 02 , 03 , and 64 are obtained from 
Zaitsev's boundary conditions; the remarkable results ob- 
tained by Eschrig Efl] arc 



61 = 



Db2 + Rbs + 046263 



03 = 



1 + 04 (Dbs 
Da4 + Rai -\ 



- i?62) ' 

010462 



1 + 62 (Dai + i?04) 



64 = 



Dai + Ra4 + 010463 
1 + 63 (L'04 -t- Rai) 
Db^ + Rb2 + O16263 
l + oi {Db2 + Rb^) ■ 



Notice that D — \ gives: oi = 02, 03 = 04, and D — Q 
gives: 01=03,02 = 04, which correspond to the com- 
pletely transparent and completely reflective cases re- 
spectively (the same argument also holds for the 6's). 
The values of the o's and 6's at other points are calcu- 
lated using Eq. (p|). 



2. Eff'ect of surface roughness 

The effects of surface roughness are accounted for by 
introducing a thin layer of impurities (elastic scatterers) 
of width d ||5C|-|5^. In the Born approximation, the Eilen- 
berger equations in the layer are written as 

d - r 

VF • ^G,^(vF,r) + w„f3 + A„(vi^,r),Gc^(vF,r) 



0; 



A. 



A, 




(Cll) 



where 



u„{x) =UJ+ 



2t 



A,ix) = A- 



2t 



T — vpl, and I is the mean free path inside the layer. 
The degree of roughness is given by the ratio p — d/l 
in the limit when d — > 0, / ^ simultaneously. For 
strong scattering, the x-independent terms in the above 
expressions can be dropped to obtain the Schopohl-Maki 
transformed equations in the form 



cos t'— = 2uja — A 
ox 

cos9^ = 2ub-A* 
ox 



A*a^; 
I-A62, 



(C12) 
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where x = x/d and 

In the above mentioned Uniit, integrating the equations 
over X from to 1, and assumi ng that fuo^gu are slow 
functions of x, we can use Eq. (C4) with /i = 1 to cal- 
culate the jump of a^^ hi, across the boundary. This ap- 
proach also works for a rough surface (we put formally 
D = 0). 




0.2 0.4 0.6 0.8 1 



FIG. 17. Superconducting gap as a function of temperature. 



APPENDIX D: ORDER PARAMETERS IN A 
D-WAVE SUPERCONDUCTOR 

As mentioned in appendix ^ the order parameter in a 
superconductor is related to the anomalous Green's func- 
tion via the self-consistency equation ( A4). In general the 
interaction potential V{vp,-v'p) in ( A4) can have compo- 
nents in different channels. Keeping only the important 
terms, we can write 

F(vF, v^,) = Vd cos 26 cos 29' + Vd- sin 20 sin 20' + Vs 

(Dl) 

where Vd, Vd', and Vs are the components of the interac- 
tion potential in the dx^-y2, dxy, and s channels respec- 
tively. (In these expressions we assume x |j a and y \\ b, 
since the orientation of the order parameter is linked to 
the local crystal axes of the system a, b.) The order pa- 
rameter in this case can also have mixed symmetry 

A(vi?,r) = Ad(r)cos20 + Ad'(r)sin26'-h A,(r) (D2) 



Substituting into (A4) we find 



A<j(r) = AnXdT 
A,(r) = 27rA,r' 



UJ>0 



dO ^ 

ZTT 

dB_ 

2tt-'^ 



Ir) cos 261 
0,r) sin 20 



UJ>0 



dO 
2^ 



Ue,v) 



(D3) 
(D4) 
(D5) 



where Xd^d' — N{Q)Vd,d' /'^ and — N{Q)Vs are dimen- 
sionless interaction constants. These equations are used 
in the self-consistent calculation of the components of the 
order parameter. 

In a homogeneous superconductor, the anomalous 
Green's function f^^ takes the simple form 



\m\' 



(D6) 



Keeping only the (i^2_j,2 component of the order param- 
eter we find 



Ad(r) = Ad4^r^ / 



dO Ad(T)cos2 26i 
2^ _^ Ad(r)2cos2 20 



(D7) 



where Uc is the cutoff frequency. In our numerical cal- 
culations we take ujc — IOtt. At T — 0, the right hand 
side of Eq. (D7) diverges as A^ — > 0. Thus, no mat- 
ter how small A^ is, t here exists a finite value for A^ 
that satisfies Eq. (D7). There also exists a finite tem- 
perature Tc, below which A^ is nonzero. This is not 
true for the subdominant order parameters, because the 
presence of the dominant order parameter, at T < Tc, 
introduces a cutoff which removes the divergence. As a 
result, even below Tc2, the bulk of the superconductor 
contains only one symmetry of the order parameter (the 
dominant one) even if the interaction potential supports 
other symmetries as well Near a (llO)-surface of 

a d-wave superconductor on the other hand, the dom- 
inant order parameter is suppressed and subdominant 
gaps may appear. The resulting mixed order parameter 
(e.g. d + id' or d + is) breaks time reversal symmetry and 
results in spontaneous currents near the surface. This 
behavior is clearly demonstrated in Fig. [l^a. 
One can show that A^ is related to Tc by 



A-i =ln- 

d Tc 



2^ry i 

^ 1,1 



(D8) 



LU>0 



The same equation holds for other symmetries of the or- 
der parameter. In our numerical calculations we use (DS) 
to find Xd' and As from Tc2 and Tcs respectively. 

At T = 0, and in the weak coupling limit A^ ^ 1, it 
follows from (D7) that 



Arf(O) = 2tJc/3e-l/^^ ln/3 :=:ln2- 1/2 w 1.21. (D9) 
The critical temperature Tc on the other hand is 



Tc = -ujcie 



ln7 = 0.577, 1.78. (DIO) 



Thus, Ad(0)/Tc = 7r/3/7 w 2.14. In terms of T^, Eq. (^) 
can be presented in the form : 



2Tr 



T ^ f dO 

UJ>0 \ n 



cos2 29 



1 



y/iv^ + AdiTfcos^ 29 
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In the limiting cases, the solution of this equation reads 
I A^(0)a-3C(3)1 f^'' 

Ad(T) 



A,(0) [1-3C(3)] [^-^^ 

/_327riy/^ J. Q _ XU/2 T<T 



For arbitrary te mp eratures < T < Tc the numerical 
solution of Eq. (D7) is shown in Fig. 
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